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ABSTRACT 

Keck spectroscopy of 5 proto-globular cluster candidates in NGC 1275 has been com- 
bined with HST WFPC2 photometry to explore the nature and origin of these objects 
and discriminate between merger and cooling flow scenarios for globular cluster forma- 
tion. 

The objects we have studied are not HII regions, but rather star clusters, yet their 
integrated spectral properties do not resemble young or intermediate age Magellanic 

Based on observations obtained at the W. M. Keck Observatory, which is operated jointly by the CaUfornia Institute of Technology 
and the University of California. 
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Cloud clusters or Milky Way open clusters. The clusters' Balmer absorption appears to 
be too strong to be consistent with any of the standard Bruzual and Chariot evolution- 
ary models at any metallicity. If the Bruzual and Chariot models are adopted, an IMF 
which is skewed to high masses provides a better fit to the proto-globular cluster can- 
didates data. A truncated IMF with a mass range of 2 — SM© reproduces the observed 
Balmer equivalent widths and colors at ~ 450 Myr. 

Formation in a continuous cooling flow appears to be ruled out since the age of the 
clusters is much larger than the cooling time, the spatial scale of the clusters is much 
smaller than the cooling flow radius, and the deduced star formation rate in the cooling 
flow favors a steep rather than a flat IMF. A merger would have to produce clusters 
only in the central few kpc, presumably from gas in the merging galaxies which was 
channeled rapidly to the center. Widespread shocks in merging galaxies cannot have 
produced these clusters. 

If these objects are confirmed to have a relatively flat, or truncated, IMF it is imclear 
whether or not they will evolve into objects we would regard as bona fide globular clusters. 

Subject headings: galaxies: elliptical, galaxies: nuclei, galaxies: individual, globular 
clusters: general 
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1. Introduction 

The discovery by Holtzman et al. (1992) of proto- 
globiilar cluster candidates in NGC 1275, the pe- 
culiar cD galaxy at the center of the Perseus clus- 
ter, was regarded by many as an important piece of 
evidence in favor of the merger model for elliptical 
galaxy formation (Schweizer 1987, Ashman & Zepf 
1992). The Ashman & Zepf model implies that glob- 
ular clusters are formed with high efficiency during a 
major merger of two gas-rich spiral galaxies. It has 
also been suggested (e.g. Fabian et al. 1984) that 
globular clusters can form efficiently in cooling flows. 
NGC 1275 has one of the largest known cooling flows 
(Fabian et al. 1981, Allen & Fabian 1997) and it has a 
Seyfert-type nucleus with strong Ha emission, as well 
as a large molecular and neutral hydrogen content. A 
number of features including its irregular shape may 
be indicative of a merger history (e.g. Zepf et al. 
1995). 

Minkowski (1957) discovered the existence of two 

velocity systems in the galaxy spectnim. The low- 
velocity (LV) system, comprising filaments of ionized 
gas whose radial velocity is similar to that of the main 
body of NGC 1275 (i.e. 5264±11 kms'S Strauss et 
al. 1992), is generally thought to be a product of the 
cooling flow. Lazareff et al. (1989) detected signifi- 
cant amounts (^ 6 x 10^ M0) of CO associated with 
the Ha filaments, perhaps revealing the source of raw 
materials (molecular gas) from which stars/clusters 
might form. The association of the molecular gas 
with the ionized filaments was inferred from the CO 
and Ha spectral line parameters (central velocity, line 
width and integrated intensity) which displayed sim- 
ilar changes between the two locations corresponding 
to the radio beam positions of the CO study. 

The CO gas was found to have such a low veloc- 
ity dispersion that it should collapse into stars be- 
fore it virializes in the gravitational potential of the 
galaxy. According to Lazareff et al. (1989), this im- 
plies cooling times of < 3 x lO'' yr which, with 6 x 
10^ Mq of material, leads to a star formation rate 
of > 2OOM0yr-i. Allen & Fabian (1997) derive an 
even larger fiow rate {456~^g^ Mq yr~^) from X-ray 
emission measurements. Far UV colors (Smith et al. 
1992) indicate that there are essentially no high mass 
stars (> 5 Mq) in NGC 1275 so a continuous cooling 
flow would have to produce stars with a steep IMF 
(skewed to low masses) or else star formation must 
have stopped at least 80 Myr ago. The 2-dimensional 



spectroscopy of Ferruit & Pccontal (1994) indicates 
that the proto-globular cluster candidates are associ- 
ated with dense ionized gas regions. For example, the 
morphology of their [Nil] isophote contour map corre- 
sponds quite closely to the distribution of the clusters. 
On the basis of their kinematics, these dense gas re- 
gions are thought to be associated with LV infalling 
filaments. There appears, then, to be evidence that 
the clusters, the LV ionized gas filaments and the CO 
gas are all associated with each other. 

There is also a system of giant HII clouds with a 
velocity of ~3000kms~^ relative to the main body 
of NGC 1275, referred to in the literature as the 
high-velocity (HV) system. The physical relation- 
ship of the HV system to the LV system is a matter 
of some uncertainty. Multr wavelength observational 
evidence (summarized by Kaisler et al. 1996) seems 
to indicate that this HV system is an infalling galaxy, 
still in front of NGC 1275. However, N0rgaard- 
Nielsen et al. (1993) argue that the major regions 
of extinction, which are probably associated with the 
HV system, are situated well within the main body 
of the galaxy. They conclude that the HV system is 
approximately halfway through the main body of the 
galaxy and/or the LV system. In any case, unless 
the proto-globular cluster candidates are very young 
(<40 Myr), it is unlikely that the HV system is associ- 
ated with their formation. Even if the HV system has 
proceeded halfway through the LV system it will have 
had only 10^ years to interact, given the high rela- 
tive velocity of the two systems (Richer et al. 1993). 

It is well known (e.g. Harris, 1991) that the specific 
frequency of globular clusters (number of clusters per 
unit galaxy light, S'at) is larger in giant ellipticals than 
in spiral galaxies and it is roughly three times larger 
for cD galaxies residing in rich clusters of galaxies 
than for "normal" giant ellipticals. A number of sug- 
gestions have been made to explain these high specific 
frequencies [see Forbes, Brodie & Grillmair (1997) for 
a summary] and many of these can be tested by ref- 
erence to NGC 1275. In particular, if globular clus- 
ter formation in cooling flows (Fabian et al. 1984) 
and/or mergers (Ashman & Zepf, 1992) are impor- 
tant factors in determining Sn, the evidence should 
be present in NGC 1275. 

West et al. (1995) were the first to suggest that the 
number of "excess" globular clusters above a "nor- 
mal" Sn correlates with the galaxy cluster's X-ray 
temperature {Tx oc Mduster)- They contend that if 
a galaxy is at, or close to, the center of a rich cluster 
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of galaxies, a high Sn globular chistcr system should 
be found. Both a central location and the cluster 
richness are necessary conditions for a high value of 
Sn- The most recent work by Blakcslcc (1996), us- 
ing a homogeneous data base comprising 5 cD and 
14 gE galaxies, indicates that Sn is indeed strongly 
correlated with cluster density as measured by veloc- 
ity dispersion, the temperature and luminosity of the 
X-ray emitting gas and, to a lesser extent, the lo- 
cal density of bright galaxies. Blakeslee's view is that 
bright central galaxies with high S n values, which are 
usually but not always cDs, are not anomalous in the 
sense that they have too many globular clusters for 
their luminosity. Rather, these galaxies are under lu- 
minous for their central position, and the number of 
globular clusters accurately reflects the local galaxy 
cluster mass density. It now appears that cluster gi- 
ant elliptical and cD galaxies cover a continuum of 
Sn values and these values depend on properties of 
the local environment. The fact that there is never 
more than one massive, high 5jv galaxy per cluster is 
consistent with the West et al. and Blakeslee results. 

NGC 1275 is located at the center of the X-ray 
emission in Perseus (Allen et al. 1992). The cen- 
tral X ray temperature and the velocity dispersion 
of the cluster are 7.5 keV (Arnaud et al. 1994) and 
1277t?|kms-i (Zabludoff et al. 1990) respectively. 
Applying these values to the relations illustrated in 
Blakeslee (1996, Figs. 5.32 and 5.27) leads to pre- 
dicted Sn values for NGC 1275 of approximately 9 
and 11 respectively, albeit with large uncertainties. 

Overall, then, a high specific frequency is to be 
expected for the old globular cluster system of NGC 
1275. The "old" globular cluster population has been 
studied by Kaisler et al. (1996) who found a Sn value 
of 4.4±1.2, much lower than expected for a central cD 
galaxy in a rich cluster (^ ll<S'7v<15) and typical of 
values found for "normal" giant ellipticals. However, 
the value they derive is heavily dependent on the val- 
ues of apparent magnitude and distance they adopt 
for this galaxy. Kaisler et al. adopted V= 11.59 (from 
RC2) and deduce a distance modulus (m— M)v = 34.9 
assuming Ho=75 kms~^ Mpc""'^ and Ay = 0.6. They 
then derived My = —23.3 and corrected this upward 
by 15% to My =—23.14 to take account of blue light 
from starburst activity. We have reworked the num- 
bers adopting V = 11.88 (from RC3) and using a 
slightly different velocity, v = 5264 km s~^ instead of 
V = 5200 kms-^ We calculate (m-M)y = 34.83 with 
the same assumptions as Kaisler et al. , i.e. Hq =75 



and Ay = 0.6. We derive My=-22.87, or My=-22.72 
with the same 15% correction (arguably an underesti- 
mate but we adopt it for consistency with the Kaisler 
et al. calculations) . Using the total number of globu- 
lar clusters inferred from Kaisler et al. (i.e. ~7900), 
this absolute V magnitude results in a specific fre- 
quency of 6.5 rather than the 4.4 quoted by Kaisler 
et al. Furthermore, the total number of globular clus- 
ters is quite poorly constrained by the Kaisler et al. 
data, which samples only a small fraction (<4%) of 
the globular cluster luminosity function. They note 
that the last data point at the faint end of their distri- 
bution (still much brighter than the peak of the lumi- 
nosity function) is surprisingly low. The total number 
of globular clusters could easily be in error by 30% or 
more, given the large extrapolation needed to make 
this estimate. A 30% increase would easily bring the 
specific frequency in line with the expected value of 
- 10. In fact, Carlson et al. (1998) report Sn ^ 10 
for the old globular clusters in NGC 1275. Hence it 
may not be necessary to invoke the formation of large 
numbers of globular clusters in the more recent past 
to produce the expected high specific frequency. 

The initial discovery of proto-globular cluster can- 
didates by Holtzman et al. (1992) from HST WFPCl 
data was later confirmed in the ground-based imaging 
study of Richer aL (1993). Holtzman ei aL claimed 
that the photometric properties of the young clusters 
are consistent with their evolution to objects resem- 
bling Milky Way globular clusters. Prom the photo- 
metric data, Holtzman et al. favored their origin in 
a merger event, whereas Richer et al. suggested for- 
mation in a cooling flow. In principle, the age spread 
of the clusters can discriminate between these difl'er- 
ent formation scenarios. Clusters formed in a recent 
merger should all be young, with ages no greater than 
perhaps a few hundred million years. Repeated merg- 
ers will, of course, confuse this simple picture. Con- 
tinuous formation in a cooling flow would lead to a 
spread in ages with some clusters possibly older than 
a billion years (Faber 1993). 

Although recent HST observations have revealed 
proto globular cluster candidates in several other 
"merger" galaxies (e.g. Whitmore et al. 1993, Whit- 
more & Schweizcr 1995, Holtzman et al. 1996), none 
have been found in a cooling flow galaxy (Holtzman 
et al. 1996). Attempts to deduce the age spread of 
NGC 1275 proto-globular cluster candidates from the 
observed color spread (Holtzman et al. 1992, Richer 
et al. 1993) have been undermined by problems with 
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the photometry (Richer et al. 1993, Faber 1993). The 
absence of a "fading vector" (with fainter clusters 
being redder) in their data was one of the reasons 
(along with a narrow color spread, later disputed by 
Richer et al. 1993) that led Holtzman et al. to sug- 
gest a narrow age spread in the proto-globular clus- 
ter candidate system. However, the simulations of 
Faber (1993) indicate that none of the probable for- 
mation scenarios is likely to produce an observable 
"fading vector". Interestingly, though, Faber's sim- 
ulations predict a significant difference between the 
merger and the cooling flow pictures in the number 
of faint clusters expected to be present at magnitudes 
just below the WFPCl detection limit. WFPC2 data 
will be relevant to this issue (Holtzman et al. 1998). 

The next step in understanding the NGC 1275 sys- 
tem is to verify that the proto-giobular cluster can- 
didates are indeed clusters rather than Hll regions. 
Once their cluster status is confirmed, determining 
their ages and radial velocities will help constrain the 
formation models. Metallicity is unlikely to be a use- 
ful diagnostic since, for very young objects (younger 
than a few Gyr) , little can be learned about metallic- 
ity from either optical colors or moderate resolution 
spectra. 

Velocity information can also provide clues about 
the origins of these clusters. It may reveal bulk flows, 
should indicate directzly whether these objects are 
associated with the high or low velocity gaseous sys- 
tems and whether the clusters are moving with the 
main body of the galaxy or with the local flow of the 
surrounding gas. 

Here we present the results of spectroscopy with 
the Keck I telescope of 8 of the proto-globular cluster 
candidates. In section 2 we discuss the observations 
and data reduction, including a detailed discussion 
of the critically important background-subtraction 
techniques. In section 3 we discuss new photomet- 
ric data from HST and compare them with previous 
estimates from HST and ground-based telescopes. In 
section 4 we measure equivalent widths in the cluster 
spectra and compare the clusters with spectroscopic 
standards. In section 5 we present the radial veloc- 
ity information. In section 6 we estimate ages and 
set constraints on the IMF using stellar evolutionary 
models and by comparison with spectra of standard 
stars and star clusters. Section 7 is a discussion of 
the implications for globular cluster formation mod- 
els. We end with a summary and conclusions. 



2. Observations and Data Reduction 

We observed 8 proto-globular cluster candidates in 
NGC 1275 in 1994 Nov/Dec, using the Low Resolu- 
tion Imaging Spectrograph (LRIS, Oke et al. 1995) 
on the Keck I 10-m telescope. The 600 1/mm grating, 
blazed at 7500 A, and a slit width of 1" provided a 
dispersion of 1.25 A/pixel and a spectral resolution of 
5-6 A. Spectral coverage was 3900-6500 A. Observa- 
tional details are given in Table 1. 

The galaxy spectrum of NGC 1275 consists of sev- 
eral components. The LV system, comprising the 
galaxy starlight and filamentary ionized gas, contains 
strong Balmer emission with broad [0111] emission 
lines at a velocity close to the systemic velocity of the 
galaxy proper (5264 km s^^). These emission lines be- 
come stronger toward the nucleus while widespread 
patchy dust appears to be significantly diminished 
very close to the nuclear region (See N0rgaard-Nielson 
et al. 1993 and Fig. |l|/Plate ??). Superimposed on 
this LV material is the HV system with strong Balmer 
emission lines that are redshifted to a velocity of 
~8000kms~^. The emission lines for both systems 
vary significantly on spatial scales that are compa- 
rable to the sizes of the clusters themselves. This 
posed a serious challenge in determining an appropri- 
ate background subtraction procedure for extracting 
the cluster spectra. Additional challenges in reduc- 
ing the data arose from instrumental limitations (dis- 
cussed below) and the proximity of the cluster candi- 
dates to the bright galaxy nucleus. All the proto- 
globular cluster candidates we observed are within 
5.3" of the galaxy center, corresponding to a distance 
of 1.9 h-i kpc at the distance of NGC 1275. An in- 
dividual cluster flux was 20% of the background flux 
for the brightest objects we observed and 4% for the 
faintest. 

Because the grating blaze angle was quite red and 
detector was relatively insensitive to blue light, the 
system throughput at the blue end of the spectral 
range was quite poor. In addition, the spectral flat 
fields, taken with internal quartz-halogen lamps, con- 
tained spatially resolved artifacts due to illumination 
and reflection off the wire-bond region of the CCD. 
These artifacts did not appear in the data. Their pres- 
ence in the spectral flats, combined with the low S/N 
in the blue, prevented flat-flelding the data on the 
blue end of the spectrum. Flat -fielding for the red end 
of the spectrum was done in the usual way. The 2- 
dimcnsional images were rectified in both the x and y 
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directions using the IRAF tasks IDENT, REIDENT, 
FITCOORDS on spectra of comparison lamps and 
then the derived transformation was applied to the 
science images using the IRAF task TRANSFORM. 
The night sky background was subtracted using the 
IRAF task BACKGROUND. What remained were 2- 
dimensional images containing the spectra of the LV 
system, the HV system and the clusters. 

The extraction procedure consisted of subtracting 
a representative galaxy spectrum from each of the 2- 
dimensional images, scaled to match the continuum 
level at the location of each cluster. Because of the 
multiple components and their variability on small 
spatial scales, the underlying galaxy light profile was 
not easily modeled from the 2-dimensional spectra 
alone. First attempts to determine the precise contin- 
uum level using only the LRIS spectral images, con- 
taining the combined galaxy, gas and cluster spectra, 
proved unsuccessful. 

Ultimately, HST WFPC2 images of the region were 
used to model the galaxy profile at the position an- 
gle of each observation. The WFPC2 images were 
rotated to match the orientation of each LRIS spec- 
tral image and then convolved with a Gaussian func- 
tion to mimic the seeing of the Keck observations. 
The WFPG2 images were then rebinned to match 
the plate-scale of the Keck spectral images. Finally, 
columns across the "slit" on the WFPC2 image were 
inspected to determine which combination of columns 
best matched the combined galaxy/gas/cluster light 
profile from the 2~dimensional Keck spectra. After 
finding the best match, the entire procedure was ex- 
actly repeated on the same WFPG2 image minus the 
clusters. The extracted light profiles from this second 
pass were used to determine the galaxy continuum 
at the location of the cluster. At most position an- 
gles the galaxy light profile revealed itself to be more 
complex than would have been possible to model us- 
ing the 2-dimensional Keck spectra alone. Not sur- 
prisingly, this technique worked better for the higher 
S/N observations than for the lower S/N ones. For 
one position angle (158.1°) we were not able to derive 
a satisfactory model of the galaxy light profile, so no 
spectra were extracted from that set of observations. 
Fortunately, the clusters observed at that position an- 
gle, H5 and 116 (identification from Holtzman et al. 
1992), were observed at least at one other position 
angle. Once the underlying continuum level was de- 
termined, the major uncertainty (as far as spectral 
features are concerned) was in selecting an accurate 



"representative" galaxy spectrum from along the slit 
to subtract from the combined object-|-galaxy spec- 
trum. Once the galaxy spectrum is removed a pro- 
file along the slit should ideally show only the light 
profiles of the clusters. In practice, there generally 
remained some residual galaxy light, attributable to 
an imperfect match between the representative galaxy 
spectrum and the true galaxy spectrum at the loca- 
tion of the cluster. Every effort was made to minimize 
this residual galaxy light, particularly in the region of 
the spectrum where the Balmer lines are located. 

With the general galaxy spectrum removed, the 
task of eliminating the emission lines remained. It 
was necessary to perform these two operations (galaxy 
subtraction and emission line subtraction) separately 
because spatial variations in the strengths and ratios 
of the emission lines (for both the high and low veloc- 
ity systems) do not correlate with spatial variations 
in the underlying galaxy spectrum/light profile. The 
technique for removal of the emission lines was sim- 
ply to extract a sample of the emission lines from 
either side of the cluster along the slit and find the 
weighted average of the two samples that minimized 
the emission lines in the cluster spectra themselves. 
Having extracted the cluster spectra, fiux calibration 
proceeded in the usual way. Gorrections for atmo- 
spheric extinction and Galactic reddening (see Sec- 
tion 3) were made using the IRAF/NOAO routines 
STANDARD, SENSFUNC and CALIBRATE. 

Of the eight objects observed, two (II12 and HI?) 
were not detected above the noise and one (113) was 
not sufficiently resolved from its nearest, and much 
brighter, neighbor (HI) for it to be extracted suc- 
cessfully. We will not consider these objects further. 
The spectra of the remaining 5 objects are shown in 
Fig. 1^. Note that we have obtained a combined spec- 
trum of H4 and H9 which are resolvable as two sepa- 
rate objects in the HST images but are not separable 
in the ground-based spectra. As explained above, 
the background contains several components each of 
which is spatially variable and the scales and am- 
plitudes of the variations are different. Because the 
background subtraction process could not perfectly 
remove each component simultaneously, the contin- 
uum shapes of these spectra are not reliable and emis- 
sion lines may have been under-subtracted (note that 
the background-subtraction process will not, in gen- 
eral, result in over-subtraction). Residual emission 
lines are particularly apparent in the lower signal-to- 
noise spectra. It is impossible to be absolutely certain 
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that these emission lines do not belong to the objects 
themselves but we believe they are artifacts, since the 
velocities of the emission-line residuals are different 
from the cluster velocities by several hundred km/s. 
Moreover, these LV emission lines can be removed by 
adjusting the background mix (at the price of intro- 
ducing other unacceptable features such as a negative 
continuum and/or absorption features from the now 
oversubtracted HV system), indicating that they be- 
long to the background emission. 

The clusters' spectra do not resemble spectra of 
HII regions. All but the lowest signal-to-noise spec- 
trum, that of H2, contain absorption features and a 
significant continuum component. The emission lines 
are demonstrably residuals from the background sub- 
traction process. In the case of H2 where the con- 
tinuum component is very weak, the presence of [NI] 
A5199 effectively rules out an HII region. [NI] A5199 
is distinctly present in the spectrum of the nuclear 
region of the galaxy and H2 has the smallest radial 
distance from the nucleus of our entire sample, sug- 
gesting that the [NI] A5199 line in the spectrum of 
H2 is a residual from the imperfect background sub- 
traction. We conclude that all 5 objects are bona fide 
clusters. 

3. Photometry 

There is some disagreement in the literature about 
the magnitudes and colors of the proto-globular clus- 
ter candidates. Richer et al. (1993) noted prob- 
lems with the WFPCl photometry of Holtzman et al. 
(1992) and Faber (1993) indicated that the CFHT 
photometry of Richer et al. (1993) must also be 
incorrect at some level. The photometric studies 
of N0rgaard-Nielsen et al. (1993) and Ferruit & 
Pecontal (1994) do not agree particularly well with 
either of the above nor with each other. We elected to 
use the WFPC2 images, taken from the HST archive, 
to derive our own B and R magnitudes for the clus- 
ters. Using circular apertures and a curve-of-growth 
type analysis we estimated total F450W and F702W 
magnitudes. These were converted into standard 
Johnson B and R magnitudes using the zeropoints 
and color transformations given in Holtzman et al. 
(1995). Our magnitudes and colors are presented in 
Table 2 along with those of Holtzman et al. (1992), 
Richer et al. (1993), N0rgaard-Nielsen et al. (1993) 
and Ferruit & Pecontal (1994). The Holtzman et al. 
V magnitudes have been corrected for the acknowl- 



edged factor of 2 error, i.e. they have been brightened 
by 0.75 mag (see Richer et al. 1993). Except for 
the H4-I-H9 clusters, the Ferruit & Pecontal R mag- 
nitudes are in excellent agreement with our WFPC2 
R values. All the entries in Table 2 have been cor- 
rected for foreground Galactic extinction correspond- 
ing to E(B-v) =0.17 (Burstein & Heiles 1984) and Ab 
= 0.71, Ay = 0.54, Ar = 0.40 and A/ = 0.26. 

N0rgaard-Nielson et al. (1993) modeled the dust 
extinction law and estimated Ay and Ai for 9 re- 
gions (designated A1-A9) in the galaxy lying between 
4.7" and 23" from the galaxy nucleus. They found 
0.20<Av<0.81 and 0.19<Ai<0.69 for these regions. 
Fig. |/Plate ?? is a WFPC2 B image of the inner 
- 15" of NGC 1275. The dust region A9, which 
appears on the WFPC2 PC image, is marked on 
Fig. |l|/Platc ?? along with the 5 clusters whose spec- 
tra are shown in Fig. y. H3 is also marked but its spec- 
trum could not be extracted because of its proximity 
to the much brighter HI. Although none of the clus- 
ters appears by eye to be positioned directly in a dust 
region as dense as A9 (Ay =0.35), variations in the 
underlying background brightness can easily disguise 
the contrast with dense dust clumps, particularly for 
regions closer to the galaxy center. Interestingly, the 
three reddest clusters for which we obtained spectra 
are clearly more affected by dust than the others. HI 
is superposed on a dust filament of width comparable 
to the size of the cluster. H6 and H4 are both located 
close to the apparent edges of dust regions. Of the 
bluer objects we observed, only H5 is in a region that 
is apparently unaffected by patchy dust. 

To quantify these impressions we created a B-R 
color image using the HST data. For each cluster we 
measured the reddest color in the region adjoining the 
cluster, and the bluest color in the galaxy at the same 
radius as the cluster. The difference between the red- 
dest and bluest color represents a rough value for the 
possible internal reddening for each cluster, assum- 
ing that the bluest color represents the true under- 
lying color. These reddening estimates are included 
in Table 2. The bluest regions of the galaxy may 
themselves be reddened, and relatively uniformly dis- 
tributed internal dust would not be detected by this 
method. Therefore the clusters could still be bluer 
than the bluest color limits we suggest. 

The cluster apparent magnitude from the WFPC2 
data can be used to derive a mass estimate. We 
used the population synthesis models of Fritze-v. Al- 
vensleben & Burkert (1995) to determine how much 
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a cluster will fade by the time it is ~ 15 Gyr old. As- 
suming a present age of ~ 400 Myr (see section 6), a 
true distance modulus of 34.23 for NGC 1275, a V- 
band M/L ratio of 2 for the clusters (if the clusters 
do indeed evolve into globular clusters they will have 
roughly this M/L ratio at ^ 15 Gyr) and Galactic ex- 
tinction of E(B-v) =0.17, we derive a mass of '--^ 2 x 10^ 
Mq for the brightest cluster, HI. This is consistent 
with the value deduced by Holtzman et al. (1992). 
In section 6 we discuss the implications of color and 
reddening for age dating and formation scenarios for 
the clusters. 

We measured the sizes of all 7 single clusters on the 
WFPC2 images to determine if the clusters are re- 
solved. At the distance of NGC 1275 a typical Milky 
Way globular cluster (effective radius Re// ~ 3 pc, 
where Re// is the half-light radius,, as defined as in 
Whitmore et al. 1993) would not be resolved, even 
with the ~ 0.05" resolution of of the Planetary Cam- 
era (PC) on HST. At the distance of NGC 1275, one 
pixel on the PC corresponds to ~17 pc so the largest 
Milky Way globular clusters Re// ^ 10 pc would 
be barely resolvable. Using the method developed 
and calibrated by Miller et al. (1997) in their study 
of NGC 7252 young clusters, we have measured the 
FWHM and the magnitude difference between 0.5 and 
3.0 pixels radius apertures and we conclude that none 
of these brightest NGC 1275 clusters are resolved. 

4. Spectral Types and Equivalent Widths 

The proto-globular cluster candidates were com- 
pared with standard star spectra covering a range of 
spectral types and luminosity classes. Ignoring the 
overall continuum shape, which we have argued is un- 
reliable because of the background subtraction pro- 
cess and which may be subject to reddening effects 
(see section 3), the best visual match is with AO to 
A3 dwarfs and class III giants. The match is illus- 
trated in Fig. ^ in which selected comparison stars 
from the stellar library of Jacoby et al. (1984) are 
presented. The proto-globular cluster candidates are 
not well-represented by class I or II giant star spectra 
nor by spectra of stars of significantly cooler or hot- 
ter spectral classification. For instance, there are no 
obvious Hel lines which rules out the hottest stars. 
There is a hint of G-band contribution to the blue 
wing of II7 which suggests an A rather than late-B 
classification. Types later than about A3 seems to 
be ruled out because of their relatively strong Ca K 



lines. To quantify these impressions we have mea- 
sured equivalent widths in various A- and B-type 
main sequence stars from the compilations of both 
Jacoby et al. (1984) and Gunn & Stryker (1983). 

It is worth noting here that there exist in the litera- 
ture and in common use a number of alternative band- 
pass definitions and methods of measuring equivalent 
widths (EWs). The method employed in the Bruzual 
& Chariot (1993) models (hereafter BC93) gives sys- 
tematically larger values of Balmer EWs than the 
Brodie & Hanes (1986) method (also adopted by Zepf 
et al. 1995). While the Brodie & Hanes bandpasses 
are well-suited to measuring features in old globular 
clusters and early-type galaxies, they are too narrow 
for measuring Balmer lines in young stellar popula- 
tions where these features are relatively strong and 
broad. The algorithm employed in the BC93 mod- 
els for measuring the Balmer line strengths is better 
suited to measuring the wide Balmer lines of young 
stellar populations. However, the choice of continuum 
level in the BC93 program is determined only by the 
flux values of the end points of the continuum band- 
passes and can, therefore, be greatly affected by noise. 
This algorithm is not appropriate for low S/N spec- 
tra. For this reason we have employed Balmer line 
bandpass definitions that have been custom-defined 
for the purpose of studying young stellar popula- 
tions. These bandpasses (referred to as LS indices 
after Linda Schroder who defined them) are provided 
in Table 3. The algorithm we use to determine the 
continuum level is based on the average flux value in 
the continuum bandpasses and so is better suited to 
low S/N spectra. The LS bandpasses have been cho- 
sen to agree well with the BC93 method when applied 
to the same high S/N spectra. 

Bruzual and Chariot (1993) models reflecting changes 
described in Chariot, Worthey & Bressan (1996), were 
circulated in 1995. These newer models (hereafter 
BC95) return EW values measured from SEDs with 
10 A/pixel resolution. This resolution is insufficient 
for accurate measurement of narrow-bandpass indices 
and, in any case, the very large widths of the Balmer 
lines in young stellar populations make the use of 
such narrow-bandpass indices inappropriate. The 
more appropriate BC95 and LS indices were mea- 
sured on stellar spectra from the libraries of Jacoby 
et al. (1984) and Gunn & Stryker (1983) to produce 
a plot of equivalent width as a function of spectral 
type (Fig. ||). The Jacoby et al. spectra were first 
smoothed to match the 10 A/pixel resolution of the 
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Gunn & Stryker spectra. As demonstrated in Fig. ^, 
the agreement between LS and BC95 indices is excel- 
lent for H/3 and H5. The small systematic offset in the 
values of H7 is due to the fact that the BC95 algo- 
rithm only uses the red side of the line when comput- 
ing this equivalent width. This suppresses any effect 
from of the G-band, which is present just blueward 
of H7. Both BC95 and LS indices are relatively in- 
sensitive to spectral resolution. 

The horizontal lines in Fig. ^ represent the LS EW 
measurements of HI, for which we obtained the high- 
est signal-to-noise data, and the arithmetic mean of 
4 clusters. H2 was excluded from the mean since 
the Balmer line emission from the galaxy was clearly 
under-subtracted in that spectrum. H/3 is subject 
to the greatest background subtraction uncertainties 
of all the Balmer lines in our spectral range because 
it is so strong in emission in the background spec- 
trum. There is ample evidence in Fig. H, especially in 
the lower signal-to-noise spectra, that H/3 emission 
from the background has been under-subtacted. The 
residual emission serves to depress the mean EW for 
H/3. Note that, in general, the background subtrac- 
tion process was executed in such a way as to result in 
either correct or under-subtraction of the background 
emission lines rather than their over-subtraction. 

Since the cluster spectra are integrated spectra, a 
perfect match with a single spectral type is not to be 
expected. Moreover, there is considerable scatter in 
Balmer EW measurements for different stars within 
a given spectral type, although these differences are 
generally comparable to our cluster measurement er- 
rors. However, the EWs confirm the visual impres- 
sions of the best matching spectral type from Fig. ^ 
and Fig. |^. Taking account of the absence of strong 
Ca K, which rules out later A-types, we conclude that 
the NGC 1275 proto-globular cluster candidates are 
dominated by early A-type stars with a best match 
at perhaps AO or Al. This conclusion is supported by 
the stellar evolutionary modeling described in section 
6. 

Table 4 gives the LS EWs of H/3, H7 and US for 
the NGC 1275 proto-globular cluster candidates. For 
completeness, in Table 5 we provide EWs for HI de- 
rived using narrow-bandpass indices, for direct com- 
parison with the results of Zepf et al. (1995). The 
definitions of the Mg6 and Fe 5270 features are from 
Faber et al. (1985), and the Balmer line bandpasses 
are those defined by Brodie & Hancs (1986). 



5. Velocities 

Velocities were determined from each unfluxed spec- 
trum using the IRAF cross-correlation routine FX- 
COR. Velocities quoted in Table 6 are the average of 
the results of cross-correlation with 2 different, high 
signal-to-noise, A-star templates with small (abso- 
lute values of <7kms~^) heliocentric velocities, ob- 
tained from the Jacoby et al. (1984) library of stel- 
lar spectra plus a synthetic A-star template made by 
combining spectra from a large number of different 
A-stars. The cross-correlation was made over the 
region blueward of H/3 to avoid bias from possible 
under-subtraction of the background emission lines. 
The Jacoby et al. spectra have a resolution of 4.5 
A so could, in principle, provide a velocity precision 
of ~30kms~^. The results from the three templates 
agree quite well, differences being a small (~10-30%) 
fraction of the velocity error. 

We exploit the presence of Balmer absorption in 
the underlying galaxy starlight to estimate the galaxy's 
"local" line-of-sight velocity at the locations of the 
clusters. At each cluster location the cluster -|-galaxy-|-gas 
spectrum (from which the cluster spectra were even- 
tually extracted) was cross-correlated against the A- 
star templates. Since the light from the clusters 
amounts to, at most, 20% of the total light producing 
these spectra, a cross-correlation that is restricted to 
the Balmer lines should give a reasonable estimate of 
the local bulk stellar velocity. The local velocity of the 
surrounding gas was estimated simply by measuring 
the line centers of the best Gaussian fits to the var- 
ious emission line profiles in the cluster-|-galaxy-f gas 
spectra. The local galaxy and gas velocities appear in 
Table 6. Excluding the composite object, H4-f9, the 
mean difference between the velocities of the clusters 
and the local stellar velocity is 47±144kms^^. The 
mean difference between the velocites of the clusters 
and the local gas velocity is 50±101kms~^. There 
is no basis for preferentially associating the clusters 
with either the gas or the stars. None of the clusters 
is associated with the high velocity infalling material. 

Zepf et al. (1995) found velocity differences for 
HI "relative to the emission line gas" of -120 and 
— 140 km s""'^ at to two different position angles. Our 
results also show a comparable blue-shift of HI with 
respect to the [OHI] emission lines, but not with re- 
spect to [01] A6300 or H/3. It is not surprising that the 
gas velocities yielded by these features are mildly dis- 
crepant. Due to the wavelength dependence of dust 
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scattering properties, measuring features at different 
wavelengths may be probing to different depths along 
the line-of-sight. At the location of H2 and H5, the 
[OIII] emission lines in the gas spectra appear to con- 
tain two velocity components with the weaker com- 
ponent located blueward of the main component. In 
these cases we measured the velocities of both com- 
ponents. The resulting velocities are listed in Table 
6. 

We find a velocity dispersion of 244±99 kms^^ for 
the clusters (excluding the composite object, H4-I-9), 
which, although very poorly constrained, seems to 
be somewhat lower than that of globular cluster sys- 
tems in dominant central galaxies (e.g. M87, Cohen 
& Rhyzhov 1997; NGC 1399, GriUmair et al. 1994, 
Kissler-Patig et al. 1998, Minniti et al. 1998). Note 
that the NGC 1275 clusters are all very close to the 
galaxy nucleus (within ~ 2 kpc) . The central stellar 
velocity dispersion of NGC 1275 has not been accu- 
rately measured but Lazarcff et al. (1989) derived a 
value of 290 kms~^ by modeling the mass distribution 
in the galaxy's potential well. Recall that Lazareff et 
al. found a velocity dispersion of '--^ 113 km s~^ for the 
CO gas in the center of NGC 1275, presumed to be as- 
sociated with the Ha filaments and the clusters. This 
value is marginally consistent with the cluster velocity 
dispersion estimates, given the size of the errors. 

The arithmetic mean velocity of the clusters, ex- 
cluding H4-f9, is 5329±121 kms~^, consistent with 
the galaxy velocity of 5264 km measured by Strauss 
et al. (1992) as well as the local galaxy velocity de- 
termined from the background spectrum. The error 
quoted is the standard error of the mean. 

6. Comparison with Star Clusters and Stellar 
Evolutionary Models 

6.1. Star Clusters 

Bica & AUoin (1986) present integrated spectra for 
star clusters covering a wide range of ages and metal- 
licities. Their plots of index EW versus metallicity 
confirm the expected lack of sensitivity to metallicity 
of line features in young cluster spectra. 

Table 7 is a list of the LS Balmer line EWs for 
the spectra of selected young clusters from the Bica 
& AUoin data base, as tabulated by Leitherer et al. 
(1996). The best match is obtained for clusters having 
ages of 200-500 Myr, though the NGC 1275 clusters 
have EWs which are larger than those for any cluster 
in the Bica & AUoin sample, irrespective of age or 



metallicity. This suggests that the NGC 1275 clusters 
are not well-represented by Magellanic Cloud young 
or intermediate age clusters, nor by Galactic globular 
or open clusters, a conclusion that is confirmed by 
visual inspection of the spectra of the Bica & AUoin 
objects (Fig. I). 

6.2. Bruzual & Chariot Models 

BC95 stellar evolutionary models at solar metallic- 
ity were used to generate plots of H/3, H(5 and H7 EWs 
as a function of age. These are shown in Fig. ^. The 
horizontal lines in these plots represent the EWs mea- 
sured for the highest signal-to-noise cluster, HI, and 
the arithmetic mean of 4 clusters. Again, H2 was 
excluded from the mean because significant residual 
emission lines from the background galaxy are present 
in its spectrum. Note that, as mentioned in section 
4, the mean EW of H/3 is likely to be somewhat de- 
pressed due to under-subtraction of the background 
H/3 component. This may be true for HI but is espe- 
cially relevant for the lower signal-to-noise spectra. 

Initially a Salpeter IMF (0.1- 125 Mq , X=1.35) 
was adopted and a delta function burst of star for- 
mation at t=0 years was assumed. The relatively low 
S/N of our spectra (with the exception of HI) pre- 
vents us from using the exact prescription of BC95 
for a comparison of the clusters' Balmer line strengths 
with the model predictions. The smoothing and re- 
binning to very low resolution (10 A/pixel) combined 
with the algorithm's "end-point only" determination 
of the continuum flux can produce quite spurious re- 
sults for low S/N spectra. For this reason we adopted 
the index measurement method outlined in Brodie 
and Huchra (1990) using the LS bandpass definitions 
decribed in section 4. Fig. || demonstrates that, for 
high S/N spectra, this method and the BC95 algo- 
rithm produce very similar results. For low S/N spec- 
tra this method allows a much more accurate deter- 
mination of the continuum flux and therefore provides 
the most robust comparison with the models. 

In Fig. H the lines representing HI, the cluster with 
by far the highest signal-to-noise spectrum, do not 
intersect the Salpeter IMF curves. It is possible to 
force an intersection for H/3 (known to be artifically 
depressed) and barely possible for H(5 (which is in 
the lowest signal-to-noise portion of the spectrum) 
by pushing both the EWs systematically lower to- 
wards the limit of their error bars. The value for H7, 
arguably the best-determined EW, is more than two 
standard deviations away from the maximum of the 
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Salpeter curve. The sample mean EWs for H7 and 
H(5 also do not intersect the Salpeter IMF curves. A 
shift towards the lower limit of its error bar does allow 
an intersection for US but not for H7. The H/3 mea- 
surement is included for completeness but is not reli- 
ably enough determined in these low signal-to-noise 
spectra to be a useful constraint on the models. 

The latest, albeit preliminary, Bruzual and Char- 
lot (1997) models (hereafter BC97) allow metallicity 
to be varied in discrete steps (corresponding to [Fe/H] 
= -1.7, -0.7, -0.4, 0.0, 0.4 and 0.7) while earlier ver- 
sions forced the assumption of solar metallicity for 
the system. These BC97 models were kindly released 
to us in preliminary form in advance of publication. 
The theoretical input spectra used in BC97 have a 
substantially lower resolution (20 A/pixel) than the 
empirical spectra used in BC95. This has caused (at 
the time of this writing) a systematic offset from the 
BC95 models, especially at young ages, in the sense 
that the BC97 models return higher values of Balmer 
line EWs. Nevertheless, the preliminary models can 
adequately demonstrate the range of EWs that can 
be produced by varying metallicity. Fig. ^ shows the 
Balmer line EWs versus age for the range of metallici- 
ties provided in the BC97 models. If the metallicity is 
unknown, II7 gives the tightest constraint on age, for 
ages <500 Myrs, owing to its small spread in EW for 
populations of different metallicity. For each Balmer 
line the maximum value of the EW occurs at solar 
metallicity for ages less than 500 Myr. Clearly, the 
solid curves in Fig. ^ (which represent solar metallic- 
ity) cannot be brought into agreement with the values 
measured for the clusters simply by allowing metallic- 
ity to vary. Except for the artifically depressed mean 
EW of 11/3, the Balmer line EWs of the clusters are 
larger than any returned by the BC models at any 
metallicity and agreement can only be achieved by 
systematic downward shifts of all the Balmer EWs. 

It is worth noting here that the BC93 models and 
later variants suffered from an error which transposed 
the column labels on the output of EW values of 
US and II7. (Bruzual, 1997). This has been corrected 
in our plots. Both Zepf et al. (1995) and Schweizer & 
Seitzer (1993) compare their measurements of Balmer 
line EWs in proto-globular cluster candidates equiv- 
alent widths from these models. Both note that the 
(BC93, solar metallicity) models are unable to pro- 
duce EWs as strong as those they observed. 

Having established that the large Balmer line EWs 
in the cluster spectra are not reproduced by vary- 



ing metallicity in the standard models, we now ex- 
plore the possibility that the IMF of these clusters 
is different from the standard IMF and specifically 
that it is skewed to high masses. All versions of the 
Bruzual and Chariot models allow for the adoption 
of Salpeter (1955) or Scab (1986) IMFs but neither 
of these is consistent with the observed Balmer line 
strengths. The input SEDs for the Bruzual and Char- 
lot models, as distributed, do not allow for variations 
in the slope of the IMF, nor can the range of allowed 
masses be specified arbitrarily. However, Bruzual and 
Chariot kindly provided us with an input SED for 
their models, generated originally for other purposes, 
which corresponds to a solar metallicity IMF trun- 
cated to include only 2 — 3Mo stars, i.e. it has a low 
mass cut-off at 2 Mq and a high mass cut-off at 3 
Mq. This essentially corresponds to a population of 
pure A stars. We have adopted this SED to simulate 
a flatter IMF. As is clear from Fig. ^, this model re- 
produces the strong Balmer line EWs of the clusters 
and these are consistent with each other in predicting 
an age of ^^425— 500 Myr for all the clusters (with the 
mean H/3 exception noted above). 

This age estimate is consistent with general ar- 
guments based on main sequence lifetimes for early 
A-type stars. Poggianti & Barbaro (1997) confirm 
that the maximum Balmer line EWs occur in AOV- 
type stars. AO star masses are estimated to be any- 
where in the range 2.2 - 3.5 M© (e.g. Allen 1974, 
Straizys & Kuriliene 1981). If higher mass stars are no 
longer present in significant numbers, a crude age con- 
straint can be determined from the duration of the H- 
burning phase of 2.2— 3.5 M0 stars. Such stars are pre- 
dicted to stay on the main sequence for ~ 0.2—0.7 Gyr, 
depending on the isochrones adopted and the metal- 
licities assumed (e.g. Schaller et al. 1992) 

Preliminary BC97 models for (B— R) as a function 
of age are shown in Fig. |^. Since the (B — R) colors 
measured from the empirical SEDs of the older models 
and the theoretical SEDs of the newest models appear 
to be consistent with each other for this broad-band 
color, we have illustrated the dependence at all avail- 
able metallicities for the standard IMF (0.1-125 Mq) 
models. Also shown is the 2 — 3 Mq model at solar 
metallicity. In Fig. |^ the color of HI, uncorrected 
for internal reddening, is represented by a horizontal 
line. A downward-pointing arrow indicates the effect 
of the estimated internal reddening, calculated from 
the HST images as described in section 3. The av- 
erage color of the 5 clusters in the sample, corrected 
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and uncorrected for internal reddening, is also shown. 

From this figure we conclude that the most prob- 
able age of HI, deduced from its color and the 2 — 
3 Mq model is in the range 425-475 Myr, depending 
on the extent of the internal reddening. The corre- 
sponding age range deduced from the average cluster 
colors is 400-475 Myr. Note that the age derived from 
the color is consistent with the age deduced from the 
Balmer line measurements. 

If a standard (Salpeter) IMF is assumed, for all but 
the lowest metallicity models, the age estimate drops 
to 400 Myr for the uncorrected color and less if any 
internal reddening is present. The age estimate drops 
to an untenable (in view of the spectral signature) 100 
Myr if the full internal reddening estimate is applied. 
We have argued that the majority of the clusters are 
likely to be suffering significant internal reddening. 
Moreover, significantly sub-solar metallicities can be 
ruled out for ages as young as those derived from the 
EW measurements. 

Taken as a whole, then, the cluster data appear 
to be more compatible with a truncated and/or flat- 
tened IMF than with a Salpeter IMF. In particular, 
age estimates from both EW and color measurements 
are mutually consistent without the need to impose 
systematic shifts of the measured values or discount 
internal reddening. 

Although the agreement of our cluster data with 
the 2—3 M0 model suggests that the IMF of the clus- 
ters may be relatively flat, i.e. skewed to high masses, 
it does not, by itself, set any constraints on the na- 
ture of the IMF above 3 Mq. After ~400 Myr, only 
A-stars will be left in the clusters regardless of the 
initial presence or absence of higher mass stars. How- 
ever, the characteristics of an original population of 
high mass ( > 3 Mq) stars, if any, can in principle be 
deduced because such stars cannot be allowed to dis- 
rupt the clusters prior to 400 Myr. In other words, the 
amount of mass lost from the clusters caused by the 
evolution of massive stars must be small enough that 
the clusters do not become unbound. Note too that, 
with a 2 Mq lower mass cut-off, these clusters would 
fade away extremely quickly (in ^ 10^ yr), assuming 
they manage to avoid disruption. Future modeling 
efforts will explore variations in the IMF slope and 
perhaps a variety of mass cut-offs in an attempt to 
reproduce the observations even more closely. 



6.3. Fritze-v. Alvensleben and Kurth Models 

Fritze-v. Alvensleben & Kurth (1997, hereafter 
FK) kindly compared our HI spectrum with their 
models. Interestingly, the FK models can reproduce 
the observed H/? and H7 EWs without the need to 
invoke a non-standard IMF, although the agreement 
with H(5 is less good. 

Table 8 gives the Balmer line EWs and colors mea- 
sured from model SEDs of various ages. FK models 
in the range 200 — 500 Myr provide a reasonable fit 
to the observed EWs with closest agreement for the 
350 Myr model. For an age of 350 Myr, the mass 
deduced for this cluster is ~ 10^ Mq. If the cluster is 
older the mass estimate will increase, if younger it will 
decrease. According to the FK models, a 400 Myr old 
cluster will fade by 3.4 magnitudes in 15 Gyr while a 
250 Myr old cluster will fade by 3.8 magnitudes (O. 
Kurth, private communication). The current abso- 
lute B magnitude of HI is -15.5, so in 15 Gyr it will 
have Mb ~-11.9, which is comparable to the bright- 
est globular clusters in M31, M87 and NGC 1399. HI 
is by far the brightest of the proto-globular cluster 
candidates but it is clear on the basis of their lumi- 
nosities, and hence their masses, that none of these 
objects would be defined as an open cluster. 

The unreddened color predicted by the 350 Myr 
FK model is (B-R)o=0.44. The B-R color of HI is 
0.56, uncorrected for internal reddening, and is 0.24 
if corrected for the value of extinction estimated in 
the immediate vicinity of the cluster for the dust lane 
in which HI appears to lie. The model continuum is 
too red compared to the observed spectrum. We ar- 
gued previously that the precise continuum shapes in 
the observed spectra are poorly defined because they 
are dependent on the background subtraction process. 
However, Fig. || and Fig. || indicate that the observed 
continuum shape is too red for the inferred spectral 
type, rather than too blue. In addition, no correction 
has been made for internal reddening. For these rea- 
sons, we infer that the true continuum shape of HI is 
bluer than that appearing in Fig. |^ and that appropri- 
ate corrections (for background subtraction inaccura- 
cies and internal reddening) would act to increase the 
discrepancy between the model and the observation. 
This discrepancy is in the sense that the model pre- 
dicts too much red light and is thus consistent with 
the suggestion that the cluster IMF is flatter than a 
Salpeter IMF. 

It is unclear why the BC95 and FK models give 
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such different EW results. Both use the Padova evo- 
lutionary tracks (Bressan et al. 1993) but there are 
subtle differences between the two, such as the details 
of the template spectra adopted. 

6.4. Metallicity 

Table 9 compares the strengths of Mg, Fe and Na 
features in the spectrum of HI with corresponding 
values from the 350 Myr FK models at various metal- 
licities. As expected for a young stellar population, 
the weak metal lines arc of little value in deducing the 
overall metallicity. However, despite the large errors, 
it is clear that Na is significantly overabundant com- 
pared to its solar value. The error on Mg is too large 
to judge its strength clearly but it too may be over- 
abundant. We note that Na and Mg overabundances 
have been seen in giant elliptical galaxies (Worthey et 
al. 1992) and in a small subset of globular clusters in 
NGC 1399 (Kissler-Patig et al. 1998). Worthey et al. 
suggest various mechanisms for producing the Mg/Fe 
enhancements seen in the gE galaxies, one of which is 
a flat IMF. 

7. Discussion 

If we assume that the clusters in our sample are 
representative of the system of proto globular clus- 
ter candidates as a whole, we can draw some inter- 
esting conclusions from our results which pertain to 
our understanding of the process of globular cluster 
formation and set some constraints on the origins of 
globular clusters. 

The mean velocity of the clusters is consistent with 
the velocity of the main body of the galaxy and their 
velocity dispersion is 244±99kms^^. This dispersion 
is on the low side for globular clusters in giant ellipti- 
cal galaxies. It has been argued that the clusters are 
associated with the LV filaments, which in turn are 
associated with the CO gas. Both the filaments and 
the CO gas are presumed to have been formed in the 
cooling flow. The cooling time for the CO gas is ~ 10 
Myr, very similar to the X-ray cooling time in these 
innermost regions (Mushotzky et al. 1981). Since we 
can rule out ages this young for the clusters, we can 
rule out their formation in a continuous cooling flow 
and possibly, by association (if the CO gas and the 
clusters are indeed spatially linked), the formation of 
the LV fllaments and the CO gas in the cooling flow. 

Ha images (Lynds 1970; Baade & Minkowski 1954) 
indicate that the extent of the filaments is 2-2.5' for 



the LV system and 1 1.5' for the HV system. Assum- 
ing a "normal" (King, de Vaucouleurs, or Hubble- 
like) density profile for the Ha material would sug- 
gest an effective radius of perhaps 20 30" and a core 
radius of even less. The scale of the filaments might 
be similar to the scale of the proto-globular cluster 
candidates. The cooling flow radius is ~5.5' (Allen 
& Fabian 1997) which is a factor of 60 larger than 
the radius of our sample and a factor of 20 larger 
than the entire sample of proto globular cluster can- 
didates. Moreover, a cooling flow scenario requires 
star formation with a steep IMF to account for the 
absence of observed high-mass stars in conjunction 
with the deduced > 4OOM0yr~^ star formation rate. 
If confirmed, the flat IMF (biased against low mass 
stars) deduced for the clusters would be further evi- 
dence against their association with the cooling flow. 
All this argues positively for a merger origin for the 
filaments and the clusters. 

The clusters show a remarkable degree of central 
concentration. Prom the work of Holtzman et al. 
(1992) and Richer et al. (1993) we know that aU 
the currently-identified proto-globular cluster candi- 
dates lie within ~8 kpc of the nucleus. Richer et 
al. (1993) noted that the brighter clusters were more 
centrally located than the fainter ones. Those in our 
study, i.e. the brighter candidates, lie within ~ 1.5 kpc 
of the nucleus. By contrast, the old globular clus- 
ters lie within about 30 kpc of the nucleus (Kaisler 
et al. 1996), so the old and young clusters are very 
differently distributed. This degree of central con- 
centration appears to be greater than that seen for 
metal-rich globular clusters in any "normal" ellipti- 
cal galaxy. 

It is not clear why formation in a merger would 
preferentially form clusters only in the central few 
kpc if star formation is triggered in chaotically dis- 
tributed shocks. However, it is to be expected that 
gas from merging galaxies will be channeled rapidly 
to the center (Barnes & Hernquist 1996 and refer- 
ences therein) where it can form stars very soon af- 
ter the onset of the merger event. It is unlikely that 
widespread shocks from merging galaxies could have 
produced these clusters both on the gromids of their 
central concentration and because their velocity dis- 
persion may be relatively low. If the clusters were 
formed further out in the galaxy and later fell into the 
center, a high velocity dispersion would be expected. 
In fact, it might be difficult to find a mechanism for 
moving the clusters nearer to the galaxy center since 
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dynamical friction is not likely to be effective at these 
masses, although they might conceivably have been 
associated with an infalling galaxy. In any case, we 
can conclude that they actually formed very close to 
the galaxy center. A cooling flow or central accumu- 
lation of merger gas is then more appealing in this 
regard because the gravitational potential definitely 
steepens rapidly near the galaxy center and objects 
formed in this potential would acquire the mean ve- 
locity of the galaxy, ft will be of great interest to 
obtain similar high-quality spectra of proto-globular 
cluster candidates in other interacting galaxies, espe- 
cially ones without cooling flows. 

If mergers are to contribute significantly to high Sjv 
globular cluster systems (although see Forbes, Brodie 
& Grillmair 1997 for arguments against that idea) and 
mergers produce very centrally concentrated globular 
clusters, some mechanism has to be invoked to redis- 
tribute the clusters to resemble the observed spatial 
distributions in normal cluster systems. For example, 
in M87 (Grillmair et al. 1986 and Lauer & Kormendy, 
1986) and in NGC 1399 (Kisslcr-Patig et al. 1997) it 
has been established that the globular cluster systems 
are more diffuse than the galaxy light. A similar argu- 
ment can be made for centrally concentrated objects 
formed from a cooling flow. 

The crossing time for objects with ~200kms~^ 
velocity dispersion with a radius of 8 kpc (the radial 
extent of whole proto-globular cluster candidate sam- 
ple) is ~ 40 Myr, or a factor of 10 less than their age. 
Their survival in the galaxy center must indicate that 
they are in virial equilibrium, unless the velocity dis- 
persion measurement is unrepresentative because of 
the very small sample of velocities. 

We have found evidence that these clusters have 
an IMF that is flat (skewed to high masses or biassed 
against low masses) compared to a standard IMF. The 
low-mass stellar component may be weak or possibly 
absent in these objects and, in order to avoid disrup- 
tion long enough to survive to -^400 Myr, there are 
restrictions on the high mass component. 

Interestingly, recent work by Padoan et al. (1997) 
provides a natural framework for such a result if 
star formation is driven by random supersonic flows. 
They predict that starburst regions should have flat- 
ter IMFs with a more massive low-mass cut-off be- 
cause of their higher mean temperatures. Based on 
core-radius measurements, Elson, Freeman & Lauer 
(1989) report evidence for a flatter IMF in their study 
of the core radii of young clusters in the Large Mag- 



ellanic Cloud. However, as already noted, the spectra 
of LMC clusters are not similar to those of the NGC 
1275 young clusters. 

8. Summary and Conclusions 

The clusters we have studied are not HII regions 
and their integrated spectral properties do not re- 
semble those of young or intermediate age Magellanic 
Cloud clusters or Milky Way open clusters. 

The clusters' Balmer line strengths appear to be 
too strong to be consistent with any of the standard 
Bruzual and Chariot evolutionary models at any age 
or metallicity. If the Bruzual and Chariot models are 
adopted, an IMF which is skewed to high masses pro- 
vides a more internally consistent fit to the spectral 
and photometric data. A truncated IMF with a mass 
range of 2 — SM© adopted to mimic a flatter IMF, 
reproduces the observed EWs at 425 — 500 Myr and 
the cluster colors predict ages of ~ 400 — 475 Myr for 
this mass range. 

The FK models are better able to reproduce the 
observed Balmer line strengths for a solar metallic- 
ity cluster with an age of 350 Myr and a standard 
(Salpeter) IMF, although the model c;ontinuum is too 
red compared to the observed spectrum. The con- 
tinuum discrepancy between the model and the ob- 
servation is in the sense that the model predicts too 
much red light and thus is also consistent with the 
suggestion that the cluster IMF is relatively flat. 

Key properties of the clusters, in addition to their 
possibly flat IMF, are their extremely centrally con- 
centrated spatial distribution, their velocity disper- 
sion of 244±99kms~^ and their age of ^450 Myr. 

Formation in a continuous cooling flow appears to 
be ruled out since the age of the clusters is much larger 
than the cooling time of the molecular gas (assuming 
that the molecular gas is associated with the cooling 
flow), the spatial scale of the clusters is much smaller 
than the cooling flow radius, and the deduced star 
formation rate in the cooling flow favors a steep rather 
than a flat IMF. 

A merger would have to produce clusters only in 
the central few kpc, presumably from gas from the 
merging galaxies which is channeled rapidly to the 
center. Widespread shocks from merging galaxies 
cannot have produced these clusters. If they were 
formed further out in the galaxy and later fell into 
the center (by some unspecified mechanism), a high 
velocity dispersion would be expected. 
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The crossing time of the clusters is only ^ 40 Myr. 
Their survival to ~ 400 Myr indicates that they must 
be in virial equilibrium if the observed velocity dis- 
persion is representative of the whole young cluster 
population. 

Regardless of their origin, these objects are not dis- 
tributed like old globular clusters in central cD galax- 
ies where the globular clusters are significantly more 
diffuse than the galaxy light. 

The clusters appear to have formed in a discrete 
event some 450 Myr ago. This event might have been 
induced by a merger which provided the fuel for a 
short-lived gas infall episode. With a small or absent 
low mass stellar component, these objects would fade 
away very rapidly (in ~ 10® yr), assuming they can 
avoid disruption in the later stages of stellar evolu- 
tion. More data and further modeling, especially to 
determine the slope and/or cut-off parameters of the 
IMF, are required to ascertain whether or not they 
will eventually evolve into objects we would regard as 
bona fide globular clusters. 
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Table 1 
Observational Data 



Cluster 


Slit 


Integration 


Dates 


ID 


PA 


Time 


Observed 


(1) 


(2) 


(3) 


(4) 


HI 


153.0 


900 


Dec 01 




46.7 


6300 


Nov 30 


H2 


113.5 


3200 


Nov 30 


H3 


46.7 


6300 


Nov 30 


H4+H9 


113.5 


3200 


Nov 30 


H5 


113.5 


3200 


Nov 30 




158.1 


1800 


Dec 01 




140.0 


900 


Dec 01 


H6 


46.7 


6300 


Nov 30 




158.1 


1800 


Dec 01 


H12 


77.4 


4500 


Dec 01 


H17 


77.4 


4500 


Dec 01 



("'^^ Cluster identification number, 
as taken from Holtzman et al. (1992) 
'•^■'Degrees from counter-clockwise from N-S 
(■^^ Exposure time in seconds 
(''^Date of Observation, 1994 
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Table 2 
Magnitudes* and Colors 



This F & P Richer Holtzman N0rgaard 





Paper 


(1994) 


et al. (1993) 


et al. (1992) 


et al. (1993) 






(1) 


(2) 




(3) 




(4) 




(5) 




Cluster ID 


B 


R 


R 


B 


V 


I 


V 


R 


V 


I 


HI 


19.00 


18.44 


18.50 


18.98 


18.65 


18.38 


18.41 


18.22 


18.36 


17.73 




0.15 


0.15 


0.16 


0.01 


0.03 


0.02 


0.02 


0.02 


0.01 


0.00 


H2 


20.18 


19.53 


19.54 


20.11 


19.89 


19.53 


19.51 


19.38 


19.95 


19.43 


crH2 


0.15 


0.15 


0.16 


0.07 


0.09 


0.09 


0.08 


0.14 


0.02 


0.02 


H3 


20.05 


19.35 


19.62 


20.34 


19.66 


19.48 


19.41 


19.23 






crH3 


0.15 


0.15 


0.35 


0.11 


0.04 


0.04 


0.08 


0.05 






H5 


20.57 


20.15 




20.71 


20.36 


20.05 


20.11 


20.01 


20.96 


20.56 




0.15 


0.15 




0.05 


0.06 


0.03 


0.04 


0.06 


0.03 


0.05 


H6 


20.76 


20.04 




20.63 


20.46 


19.91 


20.21 


20.02 


20.12 


19.45 




0.15 


0.15 




0.19 


0.08 


0.03 


0.06 


0.08 


0.02 


0.00 


H4+H9** 


19.72 


19.07 


20.04 


(19.91 


20.52 


20.70) 


(19.48 


19.30) 


(20.21 


19.74) 




0.15 


0.15 


0.47 


(0.06 


0.18 


0.22) 


(0.08 


0.11) 


(0.02 


0.02) 




B-R 


E(B-_R) 




B-V 


V-I 




V-R 




V-I 




HI 


0.56 


0.32 




0.33 


0.27 




0.19 




0.63 




fm 


0.21 






0.03 


0.04 




0.03 




0.01 




H2 


0.65 


0.18 




0.22 


0.36 




0.13 




0.39 






0.21 






0.11 


0.13 




0.16 




0.03 




H3 


0.70 


0.23 




0.68 


0.18 




0.18 










0.21 






0.12 


0.06 




0.09 








H5 


0.42 


0.03 




0.35 


0.31 




0.10 




0.40 






0.21 






0.08 


0.07 




0.07 




0.06 




H6 


0.72 


0.25 




0.17 


0.55 




0.19 




0.67 






0.21 






0.21 


0.09 




0.10 




0.02 




H4+H9 


0.65 


0.33 




(-0.60 


-0.18) 




(0.18) 




(0.47) 






0.21 






(0.19 


0.28) 




(0.14) 




(0.03) 





NOTES 

*Magnitudes have been corrected for foreground Galactic extinction (i.e. Ab = 0.71, Ay = 0.54, Aji = 0.40 and 
A/ = 0.26). The Holtzman et al. values have been corrected for the acknowledged (Richer et al. 1993) factor of 2 
error, i.e. .75 mags. 



**The objects H4 and H9 are unresolved from the ground. Ferruit & Pecontal (1994) note this and give a mag- 
nitude for H4-I-H9. Richer et al. (1993), however, list magnitudes only for H9, whereas N0rgaard-Nielsen et al. 
(1993) list magnitudes only for H4. Wc have placed those magnitudes on the line for H4-I-H9, since each is bound 
to include both objects. The entries for H4-I-H9 in the Holtzman et al. (1992) columns were derived by adding the 
objects' individual magnitudes. In this paper we have used an aperture encompassing both objects. 



1WFPC2, ^CFHT, ^CFHT, *WFPC1, ^NOT 
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Table 3 

Custom-defined LS Balmer Line Bandpasses 



Feature 
Name 


Blue 
Continuum 


Feature 


Red 
Continuum 




H/3ls 


4025.0-4060.0 
4225.0-4275.0 
4800.0-4830.0 


4070.0-4135.0 
4320.0-4370.0 
4830.0-4900.0 


4145.0 4180.0 
4375.0-4425.0 
4905.0-4930.0 




Table 4 

Balmer Line Equivalent Widths 


Cluster 
ID 


Position 
Angle 




H7LS 




HI 


46.7 


13.63^2.29 


14.90tl;0^ 


10.3511:1^ 


HI 


153.0 


13.38t^;38 


10.06tJ;|o 




HI Weighted Mean 




13.56tl-il 


IS 60+0.90 


io.8oil:r. 


H5 


113.5 




19.581^:^^ 


13.731^:^? 


H5 


140.0 






o rq+17.61 


H5 Weighted Mean 




OO or+3.81 


10.Z0_3 Q2 


11.3911^0 


H4+9 


113.5 


4.9ltt^^ 




-0.43tt^? 


H2 


113.5 




2.3Atlil (-36.431^1°)* 


H6 


46.7 


ll-16t9.?9 


12.55t|f? 


14 77+4.21 

J-*-' '-4.54 



Arithmetic Mean 

of All Except H2** • • • 13.37 ± 3.94 13.86 ± 1.56 9.13 ± 3.30 



* Reflects the presence of residual emission from imperfect subtraction of the back- 
ground emission line spectrum of the galaxy. Refer to Fig. 2. 

**H2 is excluded from the arithmetic mean because of the above mentioned pres- 
ence of residual Balmer line emission from imperfect galaxy subtraction. The errors 
quoted for the mean is the standard error of the mean. 
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Table 5 

Comparison of HI EWs with PT:BLisnED Data 





Position 












Source 


Angle 






Wbh 


Mgb 


Fe52 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


This study 


46.7 


11.2 


9.9 


10.9 


1.2 


1.9 


This study 


153.0 


8.7 


8.6 


11.9 


1.0 


0.1 


Zepf et al. 


0.0 


13.3 


10.7 


14.7 


0.6 


2.6 


Zepf et al. 


30.0 


8.3 


10.0 


11.5 


1.3 


1.9 



NOTE-AU equivalent width values given in A. 

Source of information 
^^^Counter-clockwise from N-S 

(3)-(5)Balmer line EWs, as defined in Brodie & Hanes (1986) 
(6)-(7)As defined by Burstein et al. (1984) 
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Table 6 

Velocities of Clusters, Galaxy and Gas 



Gas velocity based on emission lines 



Cluster 


Position 


Cluster 


Galaxy 


Ef3 




[GUI] 


[01] 


Average 


ID 


Angle 


Velocity 


Velocity 


486lA 


4959A 5007A 


6300A 


Gas Velocity 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


HI 


46.7 


5282± 87 


5221± 71 


5262 


5370 


5334 


5250 


5307±28 


HI 


153.0 


5111± 82 


5190± 67 


5251 


5397 


5388 


5261 


5327±39 


H5 


113.5 


5293±146 


5326± 66 


5367 


5429 


5370 


5337 


5379±18 


H5 


140.0 


5150±149 


5309± 68 


5346 


5289 


5349 


5349 


5336±15 


HI 


Wtd Mean 


5191± 60 


5205± 49 


5257 


5384 


5361 


5256 


5314±33 


H5 


Wtd Mean 


5223±102 


5318± 47 


5357 


5359 


5360 


5343 


5354±04 


H6 


46.7 


5208±110 


5251± 75 


5221 


5292 


5232 


5207 


5241±18 


H2 


113.5 


5694±207 


5343± 98 


5261 


5174 


5158 


5275 


5220±29 


H4+9 


113.5 


4553±104 


5263±103 


5244 


5232 


5223 


5243 


5239±04 



(^^Counter-clockwise from N-S 

(3) -(4) Measured by cross-correlation 

(5)-(8)]y[ga,sured by measuring line centers 

(^'The arithmetic mean of the four line-of-sight velocities measured for the gas. Error reported is the 
standard error on the mean. 



Table 7 

Balmer Line Equivalent Widths of Young Star Clusters 





Estimated 


Mctallicity 








Cluster ID 


Age (Myr) 


[Fe/H] ■ 




H7LS 


Hf-hs 


NGC 1847/2157/2214 


25 


-0.4 


3.40 


3.14 


3.51 


NGC 1866 


80 


-0.5 


7.50 


6.48 


6.96 


NGC 1831/1868 


200-500 


-0.6 


9.27 


8.34 


8.13 



Spectra, ages and metallicities were obtained from Bica & Alloin (1986) as 
tabulated in the data base of Leitherer et al. (1996). 
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Table 8 

Fritze v. Alvenzleben KuRTii Spectral Synthesis Model Results 



Age (Myrs) 


B-V 


V-R 


B- 


-R 


H/? (A) 


H7(A) 


m (A) 


100 


0.11 


0.16 




0.26 


8.17 


6.83 


5.16 


150 


0.12 


0.16 




0.28 


9.60 


8.08 


6.06 


200 


0.14 


0.15 




0.29 


10.52 


8.47 


6.61 


250 


0.17 


0.16 




0.33 


11.14 


8.77 


6.95 


300 


0.20 


0.18 




0.38 


11.33 


8.84 


7.02 


350 


0.24 


0.20 




0.44 


11.35 


8.81 


6.98 


400 


0.27 


0.21 




0.48 


11.05 


8.55 


6.67 


450 


0.30 


0.22 




0.52 


10.73 


8.24 


6.30 


500 


0.33 


0.24 




0.57 


10.37 


7.92 


6.05 


HI 






0.56 


-0.24 


11.12 


8.68 


9.21 



NOTES 

Note that the equivalent width values Usted here were measured on the 
model spectra using customized bandpass definitions that are narrower than 
the LS indices and therefore are comparable only to each other, not to the 
values listed in Tables 4 and 7. 

The color range for HI is defined by the color uncorrected for internal 
reddening and the color corrected for the maximum likely extinction for this 
cluster. See section 3. 



Table 9 

Comparison of Measured Metal Absorption Line Strengths with Model* Spectra of 350 Myr 

Old Clusters of Various Metallicities 



z 


Mg2** 


Mgl** 


Mgb 


Fe5335 


Fe5270 


Nal 


0.008 


0.043 


0.017 


0.67 


0.66 


0.54 


0.48 


0.020 


0.056 


0.023 


0.81 


0.77 


0.68 


0.68 


0.040 


0.079 


0.035 


1.11 


1.07 


1.01 


1.11 


HI Data 


0.001 


-0.016 


1.14 


0.58 


1.00 


5.38 


Avg Error 


0.026 


0.025 


0.77 


1.42 


0.76 


0.43 



*From the spectral synthesis models of Pritze-v.Alvensleben & 
Kurth (1997). 

**Mg2 and Mgl are given in magnitudes. All other line strengths 
are given in A. 
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Fig. 1.— WFPC2 B image (F450W) of the inner 15" of NGC 1275. The 5 clusters for which we obtained spectra are 
circled. H3, also circled, was observed but its spectrum could not be successfully extracted because of its proximity 
to HI. Also marked is the dust region A9 studied by N0rgaard-Nielsen et al. (1993). The very bright object near 
the bottom of the image is a Galactic F-star. The image is oriented such that North is 138° counter-clockwise 
from the vertical. East is 228° clockwise from the vertical. 
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4000 4500 5000 5500 6000 
Wavelength (A) 

(1) H/3 X4861 - LV system 

(2) [OIIll X4959 - LV system 

(3) OIII] A5007 - LV system 

(4) [NI] A5199 - LV system 



Fig. 2. — Keck I spectra of NGC 1275 proto-globular cluster candidates. The data have been smoothed with a 

box size of seven pixels and the resolution is 5.1 A. 

The emission lines, which are residuals from imperfect subtraction of the background galaxy/gas spectrum, are 
marked to indicate the source of the line. LV indicates that the lines are residuals from the low-velocity system 
(see text). HV indicates that the lines are residuals from the high-velocity infalling system (see text). 
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4000 4500 5000 5500 6000 4000 4500 5000 5500 6000 
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Fig. 3. — Spectra of solar metallicity main sequence stars and integrated spectra of young star clusters for 
comparison with HI. We conclude that the optical light from HI is dominated by stars with a spectral type 
between AO and A3. Of the integrated cluster spectra of Bica & Alloin (1986), HI is most comparable to the 
composite LMC young globular cluster NGC 1831/1868 with an of age of 200-500 Myr and [Fe/H] of -0.6, but is 
not particularly well-represented by this spectrum. 
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Fig. 4. — Balmer line equivalent widths of main sequence stars vs. spectral type, computed using two different 

algorithms. The open pentagons represent the equivalent widths computed with the same algorithm and bandpass 
definitions used in the BC95 models. The crosses represent equivalent widths computed using the method of Brodie 
& Huchra (1990) and the custom-defined LS bandpasses. In both cases the spectra were smoothed to 10 A/pixel to 
simulate the resolution of the BC95 model spectra. The two algorithms produce results that agree reasonably 
well. Shown for comparison are horizontal lines indicating the Balmer line equivalent widths (using the method of 
Brodie & Huchra and the LS bandpass definitions) of HI (solid line) as well as the arithmetic mean for 4 of the 
observed clusters (dotted line). H2 was excluded from the mean, since the Balmer line emission from the galaxy 
was clearly under-subtracted in that spectrum. The error bar for HI represents the photon error. The error bar 
for the arithmetic mean represents the standard error on the mean, i.e. the dispersion divided by y/N. 
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Fig. 5. — The time evolution of Balmer line equivalent widths for populations with different IMFs, produced by 
the BC95 models. Also shown are horizontal lines indicating the values of the Balmer line equivalent widths of HI 
(solid line) and the arithmetic mean equivalent widths for 4 of the observed clusters (dotted line). H2 was excluded 
from the mean since the Balmer line emission from the galaxy was clearly under-subtracted in that spectrum. The 
error bar for HI represents the photon error. The error bar for the arithmetic mean represents the standard error 
on the mean, i.e. the dispersion divided by y/N. 



27 




Age (Gyr) 



Fig. 6. — The time evolution of Balmer line equivalent widths for a simple stellar populations of various metallicities, 

produced by the preliminary BC97 models. For ages less than ^0.5 Gyr, H7 gives the tightest constraint on age, 
if the metallicity is unknown. This figure demonstrates the possible spread in Balmer line equivalent widths that 
can be attributed to metallicity differences. 
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Fig. 7. — (B— R) color vs. Age produced by the BC97 preliminary models. The color of HI, micorrected for internal 
reddening, is represented by a horizontal line. A downward-pointing arrow indicates the effect of the estimated 
internal reddening, calculated from the HST images as described in section 3. The possible range of the average 
color of the 5 clusters in the sample is shown as a shaded region delimited by horizontal lines representing colors 
corrected and uncorrected for internal reddening. 
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